Phase separation in the trapped spinor gases with anisotropic spin-spin interaction 



Yajiang Hao,^'^ Yunbo Zhang, ^ J. Q. Liang, ^ and Shu Chen^'B 

' Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100080, P. R. China 
^Department of Physics and Institute of Theoretical Physics, 
Shanxi University, Taiyuan 030006, P. R. China 

We investigate the effect of the anisotropic spin-spin interaction on the ground state density 
distribution of the one dimensional spin-1 bosonic gases within a modified Gross-Pitaevskii theory 
both in the weakly interaction regime and in the Tonks- Girardeau (TG) regime. We find that 
for ferromagnetic spinor gas the phase separation occurs even for weak anisotropy of the spin-spin 
interaction, which becomes more and more obvious and the component of mF =0 diminishes as 
the anisotropy increases. However, no phase separation is found for anti-ferromagnetic spinor gas 
in both regimes. 

PACS numbers: 03.75.Mn, 03.75.Hh, 67.40.Db 



Since Bose-Einstein condensates (BECs) of trapped al- 
kali atomic clouds were realized experimentally , many 
new regimes have been investigated extensively. The 
experimental progress on trapping cold atoms under a 
highly controllable way has opened the exciting oppor- 
tunities for studying strongly correlated atomic systems. 
When BECs are confined in a far-off-resonant optical trap 
regardless of their hyperfine state, the atomic spin de- 
grees of freedom are liberated and the spinor nature of 
the condensate can be manifested [2] . It stimulates enor- 
mous theoretical and experimental interests in studying a 
variety of spin-related properties, such as quantum entan- 
glement of spins, spin domains, etc Especially, the 
magnetism of the spinor gas has been studied by many 
authors [1, H, H, 0| • A great number of theories and ex- 
periments have shown that when the spinor gas realized 
in a magnetic trap is loaded into an optical trap, the spin 
domain will form after evolution for a period of time. 

Domain formation or phase separation in the multi- 
component BECs was also intensively investigated in the 
past years. The condensate mixture displays a novel 
phase in which phase separation occurs if there exists a 
strong repulsive interaction between the species P. For 
the spinor gases, the spin-dependent interaction is much 
weaker than the contact interaction, and thus the spin- 
dependent interaction almost has no effect on the total 
density distribution. Although each component displays 
a different density profile, no phase separation was found 
in the case of the isotropic spin-spin interaction^, [l0| . 
In this paper, we will show that anisotropic spin-spin 
interaction could result in the formation of the static 
spin domain in the spinor gases and we mainly focus 
on one-dimensional (ID) cold atom systems for both the 
weakly interacting regime and the strongly interacting 
TG regime. 

Recently, there has been tremendous experimental 
progress towards the realization of trapped ID cold atom 



systems [H [H, [11 Q. An array of ID quantum gas 
is obtained by tightly confining the particle motion in 
two directions to zero point oscillations [isj realized by 
means of two-dimensional optical lattice potentials. By 
loading BECs in the optical lattice or changing the trap 
intensities, and hence the atomic interaction strength, 
the atoms can be made to act cither like a condensate or 
like a TG gas (TgI Iv^ . The important parameter charac- 
terizing the different physical regimes of the ID quantum 
gas is 7 = mg/h^p, the ratio of the interaction to kinetic 
energy, where g is an effective ID interaction constant, 
m is the mass of the atom, and p is the density. 

Let us consider a repulsively interacting spin-1 Bose 
condensate trapped by a harmonic potential that does 
not depend on the atomic internal states V{r) = 
^[lo'^x^ + Lij\{y'^ -f z^)], where m is the mass of each 
boson, LOx is the trapping frequency along the x (radial) 
direction, and lOy — uj^ = lo±_ is the trapping frequency 
along the y and z (transverse) directions. Assuming the 
radial confinement {htOx) much weaker than the transver- 
sal one {huj±) leads to a ID configuration, in which the 
motion of the atoms is frozen along the transverse di- 
rections. In such a situation, the external potential that 
contributes to the atomic motion reads V{x) — ^uJ^x^ . 
In second quantization language, the Hamiltonian of our 
system may be expressed as 



"H — Hq + Hint + Hspin 
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where we are assuming that the interaction between 
bosons is described by a contact two-body potential, 
which may be described by a Dirac-delta function. Fi- 
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In Hspin, the 



describes the spin-spin interaction 
anisotropy parameter A is introduced phenomenologi- 
caUy to describe the anisotropy of spin-spin interaction, 
which may arise from the anisotropic magnetic dipole- 
dipolc interaction, whereas A — 1 corresponds to the 
isotropic model. When C2 < the system is in a ferro- 
magnetic regime, while for C2 > 0, the system is anti- 
ferromagnetic. We limit our discussion on < A < 1 
because only in this regime the phase separation may 
take place [l^. Here (x) {x) ) is the field opera- 
tor that annihilates (creates) an atom in the i-th internal 
state at location a;, i = +, 0, — denotes the atomic hyper- 
fine state \F = 1,TO_F = +1,0, —1), respectively. Summa- 
tion is assumed for repeated indices in the above Hamil- 
tonian and the pair of colons denote the normal-order 
product. Fx, Fy and F^ are the spin-1 matrices with the 
quantization axis taken along the z-axis direction. The 
atomic interaction constants are expressed through the 
effective ID interaction strength U0.2 with cq = U0+2U2 
and C2 = , where C/0,2 have the following relation 



2r 



(2) 
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(l-C(ao,2/rfj.)). 



Here ao.2 denotes the s-wave scattering lengths between 
two identical spin-1 bosons in the combined symmet- 
ric channel of total spin 0(2) when the cold atoms are 
trapped intensively in transverse direction with the trans- 
verse trapping frequency hu}± [13, [H, [13, [SH, d± = 
y/h/muj± and C « 1.4603. 

In order to deal with the weakly and strongly inter- 
acting regimes on the same footing, we work in a scheme 
of modified Gross-Pitaevskii theory [l^, [l^, [2^, [g^l in 
which the energy density e (p) is taken from the exactly 
solvable problem of a three-component Bose gas[l3|. It 
follows that the properties of a spinor gas are determined 
by the following spin-dependent energy functional 



£ = 



dx 



zJkl 



(3) 



A[{F,),^iF,),i + iFy\^{Fy),^]^,<i>u 



where p = J2iPi ^ J2i ^^^^ the energy density [23, 

MM 

f \ t \ [ Cop/2, 7 < 1 

e(p) = ^pe (7) = {^2^2^2/6^^ 7»1- 



The first line of energy functional (3) is made up of three 
contributions: the first one derives from the usual kinetic 
energy operator, the second one represents the additional 
term related to the inhomogeneity due to the external 
confinement V{x) [13], and the last one corresponds to 
the energy density in the homogeneous system. In the 
second and third lines, the contribution deriving from 
spin-spin interaction is involved; in particular, the term 
in the square brackets has the explicit form 



A (2poP- + 2p+po 



-pI 



P- - 2p+P- 



(5) 



At this point, we stress that, generally, the spin-spin in- 
teraction coupling (02) constant is much smaller than the 
s-wave interaction one (cq), i.e. C2 << cq. If we assume 
that there is no spin-spin interaction (c2 — 0) and no 
external trapping {V{x) = 0), then the system described 
by Hamiltonian (1) is integrable ^2d\; its ground-state en- 
ergy density has the same form as that of Lieb-Liniger 
problem |25j . 

In the system both the total atom number and the 
magnetization = J dx (F) = J dx [^*^^+ - ^*_^^] 
are conserved fo', '2^ . In order to obtain the ground state 
from a global minimization of £ with the constraints 
on both N and A4, we introduce separately Lagrange 
multiplier B to conserve M and the chemical poten- 
tial p, to conserve N. The ground state is then deter- 
mined by a minimization of the free-energy functional 
T = £ — pN — BAi. The dynamics of $i is governed by 
the coupled GPEs 
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with 



[H-B + C2 (p+ + Apo - p-)] + C2A<^l<i>*_, 
[if + C2(p++p-)]$o + 2c2A$+$_$5, (6) 
[H + B + C2 [p- + Apo - P+)\ + C2A$2$;, 



H = -^^ + V{x) + F{p) 



and 



F{p) = ^[pe{p)] = { 



dp 



.2t2 



Cop, 7 < 1 

'pV2™, 7 > 1 



(7) 



(8) 



We obtain the ground state of spin-1 BECs by propa- 
gating the coupled GPEs Eq. ^ in imaginary time. In 
each propagating step, the wave function is normal- 
ized to conserve the atomic number and by adjusting the 
Lagrange multiplier B the conservation of magnetization 
M. is assured. In our procedure the Crank-Nicholson 
scheme is used. We will determine the ground state for 
the ID spinor Bose gases trapped in the harmonic trap 
V (x) = ^muj'^x'^ both in the weakly interacting regime 
and in the TG regime. 
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FIG. 1; The density profile of ground state of a spin-1 Rb 
condensates for tlie + (solid line) component, (dashed line) 
component and - (dash-dot lines) component in the weakly 
interaction regime with m=0.2. (a) A=1.0; (b) A=0.9; (c) 
A=0.8; (d) A=0.5; (e) A=0.2; (f) A=0.0. In this figure the 
length is in the unit of a = 1.2/im. 
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FIG. 2; The density profile of ground state of a spin-1 Rb 
condensates for the + (solid line) component, (dashed line) 
component and - (dash-dot lines) component in the Tonks 
regime with m=0.2. (a) A=1.0; (b) A=0.9; (c) A=0.8; (d) 
A=0.5; (e) A=0.2; (f) A=0.0. In this figure the length is in 
the unit of a = 8.5/im. 



To investigate the effect of of anisotropy parameter 
A, we evaluate the density profile of the ground state 
of ID spinor gases for the ^^Rb (ferromagnetic) with 
flo — 102aB and 02 = lOOas {clb is the Bohr radius) 
in the harmonic trap for different anisotropy parameters. 
By properly tuning the parameters, the system may be ei- 
ther in the weakly interacting regime or in the TG regime. 
Let us first consider the specific system with the typical 
parameters of the trap lOx = 0.5kHz, lo±^ = 50kHz and 
the atomic number N = 2000, in which case the effective 
interaction strength 7 ~ 0.008 indicating that the sys- 



0.06 
0.04 

0.02 
0.00 



^ (a) 




1 1 I I 


V 


























•' V ^ 




1 1 1 1 


1 1 1 1 





4 8 12 16 4 8 12 16 4 8 12 16 
XXX 



FIG. 3; The density profile of ground state of a spin-1 Rb 
condensates for the -I- (solid line) component, (dashed line) 
component and - (dash-dot lines) component in the weakly 
interaction regime with A=0.2. (a) m=0.0; (b) m=0.5; (c) 
m=0.8. In this figure the length is in the unit of a = 1.2/im. 



tem is in the weakly interacting regime. Fig. 1 displays 
the density profiles in unites of N/a with a — ^jTi/raujx 
for different anisotropy parameter A in the weakly inter- 
action regime with m= ^ = 0.2. When the spin-spin 
interaction is isotropic, i.e. A = 1.0, the three differ- 
ent components superpose each other and they have sim- 
ilar distributions. In the case of weak anisotropy, for 
instance A = 0.9 here, the distribution has changed ex- 
plicitly. Therefore the ground state configurations have 
positive magnetization in some region but negative in 
another. Also the component (dashed line) diminishes 
as the anisotropy increases. As the anisotropy becomes 
more and more clear, the components tend to separate 
and coincide only at the boundary between them. The 
components always try to avoid each other. It is shown 
that when A — 0.2 the component disappears com- 
pletely and the density profiles of -I- component (solid 
line) and — component (dash-dot lines) exhibit in the 
form of phase separation. The corresponding density pro- 
files in the TG regime are plotted in Fig. 2 with m= 0.2. 
In this regime the parameters are tuned to ojx = 10 Hz, 
ijj±^ = 500kHz, the atomic number = 50 and the effec- 
tive interaction strength 7 ~ 15. In this case, the similar 
density distributions occur. Comparing the Fig. 1(a) and 
the Fig. 2(a) we see that the density profiles in the TG 
regime behave like that of Fermions. According to Fig. 
2(b), with the very weak anisotropy, obvious phase sepa- 
ration has occurred in the Tonks regime. Fig. 3 displays 
the density profiles for different magnetization but with 
the same anisotropy parameter, which indicates that the 
magnetization only infiuences on the ratio of the atomic 
numbers between the -I- component and — component. It 
turns out that although the spin-spin interaction is much 
small, its property affects greatly the density distribution 
of each component of the spinor gas and the phase sep- 
aration occurs more easily in the TG regime than in the 
weakly interacting regime. The ID spinor gas in the TG 
regime might provide us a good platform to investigate 
the magnetism of the cold atoms. 

For further studying the effect of anisotropy on the 
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FIG. 4: The density profile of ground state of a spin-1 Na 
condensates for tlie + (solid line) component, (dashed line) 
component and - (dash-dot lines) component in the weakly 
interaction regime with m=0.2. (a) A=0.5; (b) A=1.0. In 
this figure the length is in the unit of a = 7.4/^m. 

anti-ferromagnetic system, we consider a condensate of 
^^Na in the weakly interacting regime with oq = SOas 
and a2 = 55. las . The trap frequencies are chosen 
as ujx = 50Hz, uj± = lOkHz and N = 1000 so that 
7 ~ 0.001. The density profiles are shown in Fig. 4 
with m — 0.2 for anisotropic and isotropic case. It is 
shown that no phase separation occurs even for very 
large anisotropy and the same distribution as the case 
of isotropic spin-spin interaction displays JKi]. From the 
above results, it is obvious that the density distributions 
of every component strongly depends on its' anisotropy 
and the ferromagnetic or the anti-ferromagnetic proper- 
ties of the relatively weak spin interactions, whereas the 
total density is almost not affected by the weak spin in- 
teractions. 

Finally, we discuss the possible experimental realiza- 
tion of anisotropic spin interaction in spinor gases. It 
is well known that magnetic dipolar interactions are 
anisotropic despite the fact that dipolar interactions are 
rather weak comparing to the spin interactions. The rel- 
ative strength of the dipolar and the spin exchange in- 
teractions is estimated to be 10~^ for *^Rb and 10"'^ for 
^^Na [29]. Since the magnetic dipole-dipole interaction is 
irrelevant to the s-wave scattering length, we may tune 



the s-wave scattering length experimentally by the Fesh- 
bach resonance so that the strength of isotropic spin-spin 
interaction C2 is comparable to that of dipole-dipole in- 
teraction, and thus the anisotropic interaction becomes 
obvious. Although a condensate of strongly anisotropic 
spinor gases remains to be realized, experimentalists in 
Ref. [lOl have already successfully demonstrated the effect 
of the dipole-dipole interaction. With the present rapid 
development in the experimental manipulation of cold 
atoms, the goal of making a condensate of spinor gases 
with anisotropic spin interaction does not seem to be far- 
fetched. It is worth indicating that our approach can be 
directly applied to deal with the three-dimensional (3D) 
problem for which the mean-field theory corresponds to 
our weakly interacting theory. It follows that a 3D spinor 
gas with anisotropic spin interactions also displays phase 
separation. However, the result for the strongly interact- 
ing regime can not be extended to the higher-dimensional 
case in which no TG gas can be realized. Our work is 
helpful to understand the properties of the spinor con- 
densates and deepen our understanding of formation of 
spin domains in spinor gases. 

In summary, we have studied the density profile of 
the ground states of ID spin-1 Bose gases for different 
anisotropy parameter A. The distributions of the fer- 
romagnetic spinor gas are affected tremendously by A 
although the C2 term is very small compared with cq 
term in Eq. ([T]). Even if the anisotropy is weak, the 
distributions show obvious difference from that of the 
isotropic case. In the large anisotropy the component of 
mp = disappears and obvious phase separation occurs 
both in the weakly interaction regime and in the Tonks 
regime. And the effect of anisotropy in the TG regime 
can display more obviously with weaker anisotropy than 
the former case. However, when the spinor gas is anti- 
ferromagnetic, the distribution is no longer being affected 
by the anisotropy parameter. 
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